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By assuming the existence of a pseudo-potential smooth enough to do Hartree-Fock variations 
and good enough to describe nuclear structure, we construct mass formulae that rely on general 
scaling arguments and on a schematic reading of shell model calculations. Fits to 1751 known 
binding energies for N,Z> 8 lead to RMS errors of 614 keV with 14 parameters and 388 keV with 
28 parameters. The latter is easily reduced to a 20 parameter form at 423 keV. 
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", Calculations of nuclear masses reflect the traditional cleavage between local methods of shell model origin and global 
ones, in which semiclassical arguments play an important role. Out of the ten contributions to the 1986-1987 mass 
^\ , predictions H, two are simply numerical, three follow Garvey-Kelson lines, and five incorporate a Bethe-Weizsacker 
■ liquid drop (LD) formula as basic ingredient. The earlier work of Liran and Zeldes || is a shell model approach and 
the recent Thomas-Fermi (ETFSI) calculations of the Montreal-Brussels group || are definitely global but belong to 
Oh. a special category in that the number of parameters is truly small (9) and the LD form comes as output. 

' At present, only the local formulae of Masson and Janecke 0] can go below RMS errors of 400 keV but they need 
hundreds of parameters and cannot reach safely the drip-lines as demanded by calculations of r-processes nucleosyn- 
thesis. The droplet (FRDM) 0^],and ETFSI M mass tables are designed to extrapolate efficiently but they differ in 
the way they do it, and with RMS errors of some 700 keV, are not sufficiently precise for the more terrestrial uses. 
t-H ■ To obtain some improvements we propose to abandon the local-global dichotomy and go back to fundamentals. 
The only assumption we shall make is that there exits a nucleon-nucleon (NN) pseudopotential smooth enough to 
do Hartree-Fock (HF) calculations and good enough to explain nuclear structure. The pseudopotential cannot be 
identical to - although we expect it to resemble - the bare NN interaction. A formal consequence of the assumption 
is that we can separate the Hamiltonian as H = H m + Hm ■ The monopole contribution H m — Hf n + H' m contains a 
diagonal part 

H d m = £> w m fc (m, - 5 kl ) + b u (t k .ti - ~^5 kl ) (1) 

k,l 

- written in terms of number(m) and isospin operators(i) for the orbits k, I- whose expectation value for any state is 
the average energy of the configuration to which it belongs. (A configuration is a set of states with fixed m and t for 
each orbit). The non diagonal term H' m is such that H m is invariant under unitary transformations. Therefore HF 
variation of H amounts to HF variation of H m for closed shell and single particle (or hole) states built on them, since 
for this set (which we call cs ± 1) each configuration contains a single member. 

The multipolc Hamiltonian Hm contains all other terms (pairing, quadrupole, etc). It is therefore responsible for 
5_j ■ correlations, while H m is in charge of saturation properties. 

The H m + Hm separation is a rigorous result whose proof is sketched in |5| and given in general in 0. It was 
used in to demonstrate that once H m is treated phenomenologically, a parameter-free Hm derived from realistic 
interactions is sufficient to provide high quality shell model spectroscopy in all regions where exact diagonalizations 
are feasible (RMS errors below 300 keV in the p and sd shells, below 200 keV in the N = 50,82 isotones and Z = 
28,50 isotopes). 

In shell model calculations, the full H is not used directly but replaced by a renormalized version H = H m + 
Hm adapted to finite spaces in regions bounded by magic closures. Then Hm is taken to reproduce exactly the 
corresponding cs ± 1 states. 

Although our ultimate aim would be to discover Hm, in this paper we propose a simpler, and probably necessary 
first step which is the determination of the gross features of H m by including it as basic ingredient in a mass formula. 
Let us assume then that we have some H = H m + Hm , ready for shell model calculations : 

H = H A kim k mi + B u t k .ti + H M (2) 

which is taken to be a strictly two body force by including in it the kinetic energy after elimination of the center 
of mass contribution. However, for simplicity we have omitted in eq.(2) the counter terms in Ski in eq.(l) (we have 
checked that their effect is negligible in the fits). 

To obtain some clues about H m we reduce it to separable form by diagonalising the A^i and Bki matrices. 
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For Aki we have 



^A kl m k mi =^2e^(^2 k m k fk fl ) 2 (3) 

k,l fJ- 

and borrow from Q the result that a realistic force produces a strongly dominant eo, whose eigenvector M has 
amplitudes 

f ka = [{p+l){p + 2)}- 1/2 =D p - 1 ' 2 (4) 

where k belongs to the p-th harmonic oscillator shell of degeneracy D p . 

Let us study this term. Setting m p — n p + z p , where n, z are number operators for neutrons(z^) and protons (it), 
filling shells p v and p w up to some p /„ and p / w Fermi level, we find : 

< M >= £ i - i[(37V) 2 / 3 + (3Z) 2 / 3 ] (5) 



/Dp 
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where we have approximated \/ (p + l)(p + 2) w p + 3/2, and used TV = ^ n Pi/ = ^(p„ + l)(p„ + 2) = |(p/„ + 2) 
and Z = |(p/„ +2) 3 . 

The eigenvalue eo must behave as a typical two body matrix element, which for a realistic force goes as 

V(u))kl mn — V(u)o)klmn + 0(uj ) (6) 

a result from ref. |Q,but adding an 0(u> 2 ) correction warranted for large oscillator constant u>. Then we know that 

huj = s 4o^-V3 + (A- 2 / 3 ) (7) 

< r 2 > 

where the leading term is the classical result Q obtained with a standard mean square radius < r 2 >= 0.86A 2 / 3 . The 
0{A~ 2 / Z ) correction comes because the light nuclei are larger than the standard estimate. By combining eqs.(6) and 
(7) we obtain the scaling law T(A) = TA^ 1 / 3 + jA~ 2 / 3 , or its more flexible generalization 

T(A) = (r/R)(l-p(T)/R)) , R = Rl/AV* (8) 

that we shall use for all amplitudes T(A) affecting operators T. The form of R is left free to allow for the better 
estimate: < r 2 >= 0.90i? 2 , with R c = A 1 / 3 (l - £ (2T/A) 2 ) 1 / 3 , £, = 0.42 §. 
For the leading monopole term we have from eqs.(5) and (8) 

TT(A) = C£™ P /VD;) 2 eo = 0(A) + 0(A 2 ' 3 ). 

This remarkable object goes asymptotically as volume plus surface LD terms and at the same time produces strong 
magicity at the harmonic oscillator (HO) closures, as can be checked by plotting T (or more precisely its expectation 
value ). To obtain the usually observed extruder-intruder (EI) closures we have to add spin-orbit (SO) effects and we 
rely on the diagonal construction (3) to propose a term orthogonal to M of the form S = ^2, SpNp -1 , where N p is a 
normalization to be determined and 

S p = pmj P — 2m rp = (l.s) p . (9) 

Here (refer to Fig.l) jp is the largest orbit in the p-th shell and rp regroups all the others. (l.s) p is the operator 
that produces the same splittings as (l.s) p and then collapses the r-orbits to their centroid value. The rationale for 
considering only two types of orbits is clear from Fig.l : we want to give top priority to shell formation (i.e. the cs 
part of the cs + 1 set ).The combinations of m k operators other than nij and m r will contribute to subshell effects 
that we incorporate in Hm and treat later. 
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FIG. 1. HO and EI major shells. 

Since M and S are symmetric combinations of (properly scaled !) m p and S p operators, the only other contributions 
we can include consistently must be symmetric in p, i.e., sums of m 2 and s p , again properly scaled. The arguments 
are exactly the same for the t operators. 

Calling Tfip — 7ip~\~Zp : tp — \tl p %p\i Sp = p{rijp-\- Zjp 2(rtrp~\~ Zrp) , and St p — p\iijp Zjp\ 2\n r p z r p\ 



we introduce the variables 



MA p = ^= , SA P = Sp , (10a) 



D p ' P 2(p+l) ' 
MT p = -^= , 5T p = -%- . (10b) 



fD- p ' P 2(p+l) 

In general A = N + Z and T = \N — Z\/2, but in combinations as above, they mean "isoscalar" and "isovec- 
tor" respectively, while M stands for "master" and S for "spin-orbit". The first part of Table 1 gives the 12 
possible symmetric quadratics obtained with operators (10) . F here stands for "full" and P for "partial". All 
these operators will be affected by coefficients that scale as T(A) in eq.(8). The FCT term is chosen to go as 
an ordinary (l.s term with £ = 0(A~ 2 / 3 ) [0. The other scalings then follow by symmetry. It could be ar- 
gued that the S operators in eq.(10) should carry an extra 0(p -1 / 2 ) factor. This uncertainty is of little con- 
sequence. On the contrary the factor is important and we shall let the fits decide in favour of a = 1/2. 



TABLE I. The operators f (called T here) in H m , H a and H d 



H m (T=\N- Z\/2) 
FMA = (J2 MApf 
PMA = J2(MA P ) 2 D^ 
FAIT = (J^MTp) 2 
PMT = Y.iMTpfDp' 
4T(T+ 1)A~ 2/3 



H a (n — D v — n, z = — z) 
SQ+ = 2(nn) 2 D 2 f~ 3 + 2{zz) 2 Df' 3 



H d (n'=n- JU, ft' =n + JU), 
QQ l [ = {n'n'D- 3/2 ± z'z'D^ 3/2 ) 2 



a = 1/2 

FSA = (£ SA p f 

PSA = J2(SA P ) 2 D^ 

FST = (£ ST V ) 2 

PST = *£(ST P ) 2 D% 

V p = -mod(N, 2) - mod(Z, 2) 

52 = nnD' 1 + zzD' 1 
SQ- = ^nnzziDvD^) 13 -'^ 2 

(z' = z- JU, z' = z + JU) 

QQ± = QQ0±QQ1 



R c = [A(l-£(2T/A) 2 )] 1/3 C = 0.42 

FCA = J2 MA p Dp 1/2 SA p Dp 1/2 

PC A = J2(MAp)(SAp)D p '- 1 

FCT = J2 MTpDp 1 ' 2 J2 ST p Dp 1/2 

PCT = J2( MT p)(STp)D p '- 1 

V c = [~Z{Z - 1) + .76[Z{Z - l)] 2/3 ]/i? c 

S3 = nn(n - n)Dl^~ l) + zz(z - z)Dl^- 1] 
= 1/2 

JU = 4 
DK = 16. 



In addition to these terms H m , includes standard pairing (Vp) and Coulomb (V c ) contributions as well as a 4T(T+ 1) 
term whose presence is necessary 

If H m has ensured shell formation, mostly of EI type, the variables that become important in modelling configuration 
mixing are n v and z v , the number of valence particles in EI spaces of degeneracy D v and D T (see Fig.l). If we assume 
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a-priori the right closures and the right boundaries between spherical and deformed nuclei, we now that extremely 
good fits can be obtained with monopole-like forms B . 

Let us see how in the present formulation Hm is expected to provide naturally good boundaries, once H m has 
provided good closures. The second part of Table 1 lists the four possible operators whose appearance is guaranteed 
under very general circumstances for spherical nuclei. To a large extent they account for subshcll effects. They are 
discussed in detail in g. The D@ factors reflect scaling uncertainties to be resolved by the fit. 

Nilsson diagrams indicate that the onset of deformation is associated with the interruption of normal spherical 
filling by the promotion of 4 neutrons and 4 protons to configurations in the next HO shell (i.e. including the j orbit 
as well as others) ||. Calling qi the quadrupole moment of these intruders and q that of the orbits left behind, we 
expect a gain in energy of the form 

(q + q,) 2 =q 2 + q 2 + 2qq t 

. In the third part of Table 1 we give the expressions for QQO and QQ1 representing the two possible q 2 contributions. 
Since qi is a constant, its effects are included in the DK term mainly meant to correct the estimate of monopole loss 
coming from H m . The qqi terms are conceptually important since qi provides the effective quadrupole strength that 
will drive the lower orbits. However, we have left them out of the table because the information and the few keV 
they bring do not justify bothering with 4 extra parameters. Conversely JU is introduced explicitly to stress that by 
varying it, JU — 4 turns out to be optimal. The form of q is equivalent to equidistant Nilsson orbits and it is scaled 
so that deformation energies have a standard A 1 / 3 behaviour. 

All operators (except V c ) are affected by factors of type (8) with £ = 0.42 in R c . Energies (taken to be positive) 
are given by the expectation values 

E(N, Z) =< H m > + <H s >(l- S d )+ <H d >S d 



= max(< H m > + < H s >, < H m > + < H d >) (11) 

the lowest possible orbits are filled for spherical nuclei (Sd — 0), while for deformed ones (Sd = 1), JU particles 
are promoted to orbits j. The calculations are conducted 



TABLE II. Parameters of the Up, 28p, 28p* and 20p fits. (V p and V c given in Table 3). 28p* fit uses R = A 1/a . 



f 


FM+ 


PM+ 


4T(T+1) 


FS+ 


FC+ 


S3 


DK 


QQ- 


PS+ 


PS- 


FS- 


FC- 


SQ- 


QQ4 


ri4 


9.33 


-0.602 


-36.08 


0.44 


3.27 


0.45 


-10.0 


6.53 














p(T)14 


0.78 




1.40 






4.76 


4.41 
















T28 


9.51 


-0.79 


-36.51 


5.19 


-15.15 


0.56 


-36.1 


21.6 


-0.7 


-0.1 


1.0 


-30.01 


0.4 


3.6 


P (r)28 


0.75 




1.40 


4.24 


3.71 


4.73 


4.75 


3.54 


5.34 


4.47 


3.61 


3.51 


4.73 




T28* 


9.55 


-0.82 


-31.83 


5.26 


-26.75 


0.57 


-36.9 


19.0 


-0.7 


-0.1 


1.05 


-26.75 


0.4 


2.6 


P (r)28* 


0.75 




1.33 


4.09 


3.55 


4.77 


4.71 


2.72 


5.21 


4.24 


3.38 


3.59 


4.79 




F20 


9.67 


-0.95 


-34.77 


5.95 


-20.73 


0.54 


-44.2 


27.5 


-0.7 


-0.16 


1.07 


-20.73 


0.4 




P (r)20 


0.825 


.825 


1.35 


3.9 


3.9 


4.65 


4.65 


2.87 


4.65 


4.65 


3.9 


3.9 


4.65 





by initializing Sd, fitting E(N, Z) in the first equality of eq.(ll) to the 1751 mass values for N, Z > 8 in the latest 
compilation |lC| ], then resetting Sd through the second part of eq.(ll) and iterating until convergence. Table 2 con- 
tains results for 3 fits with 14, 28 and 20 parameters (14p, 28p, 20p), whose RMS errors are 614, 388 and 423 KeV 
respectively. The notation for the operators is XX± = XXAzkXXT and we have preset JU = 4, a = /3 = 1/2. The 
14p fit is special in that adding one parameter gains at most 50 keV, while substracting one costs 100 keV for p(DK), 
140 keV for either FS+ or FC+ and more than 300 keV for any other choice. 

In 28p, we have added two QQ+ and two SQ— parameters for a gain of 60 keV, and 10 parameters associated with 
spin orbit effects that bring in 150 keV. With 28 parameters, Table 2 and 3 compare the £ = 0.42 (28p) and £ = 
(28p*) choices for R c . 

The 20p fit is a variation of 28p in which QQ+ is excluded and five groups of operators: (FM+,PM+), 
(4T(T + 1)), (QQ-), (FS±,FC±), (PS±,DK,SQ~) are constrained to have a single p(T) per group. 
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TABLE III. Asymptotic forms of the fits compared with a pure LD form (6p). (T2 is for 4T(T+1)) 





A 


-A 2 ' 3 


-T2/A 


T2/A i/3 


v P 


Vo 


RMS 


6p 


15.49 


17.79 


28.64 


40.23 


5.22 


.705 


2.544 


Up 


14.95 


12.43 


26.73 


36.21 


5.17 


.696 


0.614 


28p 


14.95 


12.14 


27.26 


37.51 


5.21 


.699 


0.388 


28p* 


14.98 


12.21 


28.50 


40.66 


5.25 


.699 


0.409 


20p 


14.99 


12.36 


27.24 


37.30 


5.21 


.699 


0.423 



Subshell effects. The H s operators are largely devoted to mock the energy patterns generated by subshell structure 
The S2 operator is easily absorbed in H m and the heaviest task goes to S3. Parametrization (8) is now a convenient 
tool unrelated to its original derivation (eqs.(6) and (7)). This can be detected by an anomalously large p value 
leading to a change in sign of S3 at R = p, i.e. A > 100, the region in which j orbits can start filling before r orbits 
are full. In general, a large p(T) indicates that the T operator is adding to its specific job some subshell corrections 
(e.g. p(QQ—) in 28p, for a gain of 25 keV, and p{DK) (= p(S3) !) for a gain of 100 keV already mentioned). Given 
that SO is the very origin of subshell structure it comes as no surprise that S-type operators are easily contamined 
by large p ratios. 

LD parameters and radii. Table 3 compares the parameters of 6p (pure LD) with those obtained by expanding 
eq.(5) and similar ones for FM and PM+: the combinations 1.717[(FM+) + (PM+)] and 
4T(T + 1) + (0.382 + 1.145£)[(FM+) + (PM+)] 

become the asymptotic coefficients of A and 4T 2 respectively. The factor in £ comes from the R denominators. It 
doubles for the 4T 2 /A 4 / 3 term. The gain brought about by the use of a more precise form of R c is no doubt significant 
as shown by the case of 28p. 

It should be noted that with a = 1/2 scaling, PM+ is a volume term, while for a — it is pure surface and 
p(FM+) becomes very small, a disturbing result that justifies the a = 1/2 choice. 

It is remarkable, also, that the radius extracted for V c in Table 3, r m 1.235i? c , is very close to the fitted r w 1.225i? c 

B- 

The FMT - 4T(T + 1) puzzle. The presence of the 4T(T + 1) term is demanded by the fits below 650 keV, while 
the omission of FMT is possible down to the 470 keV level. The puzzle is that it serves no purpose to treat FMT as 
a free parameter: the efficient combination is FM+ — FMA + FMT. Renormalisation effects and the cancellation 
of kinetic and potential energies are the only guesses we can propose for the emergence of T(T + 1) and FM+ as 
leading operators. 

Deformed nuclei. To avoid unphysical results we allow 5d = 1 in eq.(ll) only when DK > 0, i.e. for A > 100. The 
number of S = 1 cases is (423,367,329) for the (Up, 28p, 20p) fits, and the RMS errors of (594,302,385) are smaller 
than those for the spherical nuclei in each case. The experimental ground state bands of the 5d — 1 nuclei show 
rotational features. 

To conclude: by replacing the global-local alternative by the monopole-multipole separation of TL we can find good 
mass formulae. Better ones will come once we learn more about 7i m - 
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Pearson, A. Poves, P. Quentin, A. Sobicewsky, F.K. Thielemann, F. Tondeur and N. Zeldes for useful exchanges. 
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